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ABSTRACT 
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the  high  fluid  velocities.  Certain  geometric  design  variables  of  the  gears  were 
seen  to  affect  the  peak  velocities  for  each  case.  The  variables  most  significantly 
affecting  the  peak  velocity  appear  to  be  the  drive  ratio  and  the  face  width.  The 
high  velocities  may  contribute  to  the  noise  generated  during  meshing  of  gear 
teeth  due  to  the  jet  noise  as  a  result  of  the  high  velocity  jets  impinging  on  the 
enclosures  surrounding  the  gears  and  the  formation  of  shock  waves  at  the  exit 
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Chapter  1 


INTRODUCTION 


Smooth,  quiet  running  gears  axe  desirable  in  many  applications.  Much  has 
been  done  to  identify  possible  sources  of  noise  in  gears.  It  is  recognized  that  the 
gear  meshing  action  is  one  of  the  most  important  sources  of  noise  in  high  speed 
machinery.  Meshing  noise  is  generated  by  the  nonuniform  transfer  of  torque 
between  gears  due  to  geometric  defects  in  tooth  profiles,  gear  shaft  vibrations, 
and  bending  of  the  teeth  under  load  [l].  A  considerable  amount  of  work  has 
been  done  and  is  currently  in  progress  to  predict  and  lower  the  effects  of  these 
sources.  Ishida  and  Matsuda  [2]  studied  the  effect  of  friction  noise  on  gears 
by  examining  pitch  circle  impulse  noise  simulated  by  two  disks.  The  same 
authors  also  studied  gear  noise  due  to  surface  roughness  [3].  Fukama  et  al.[4] 
analyzed  noise  mechanisms  of  axial  and  radial  vibration  in  spur  gears.  Badgley 

[5]  computed  the  flow  of  vibration  energy  within  a  gear  drive  train  as  a  method  to 
understand  the  noise  problem  and  the  technology  needed  to  alleviate  it.  Laskin 

[6]  showed  a  procedure  to  predict  gear  noise  by  compiiting  the  mechcinical  energy 
from  the  sinusoidal  displacement  and  dynamic  tooth  force  and  then  determined 
the  fraction  of  energy  radiated  in  the  form  of  acoustic  energy.  These  studies  have 
focused  primarily  on  the  physical  contact  between  the  meshing  gears  as  a  source 
of  noise.  One  possible  source  of  noise,  which  has  not  been  thoroughly  analyzed, 

is  the  flow  of  air  and  lubricant  from  between  meshing  gear  teeth.  During  a  high 

^ _ 

speed  mesh,  the  time  duration  of  the  arc  of  approach  and  recession  can  take  on  the 
order  of  100  microseconds.  During  this  cycle,  the  air  and  lubricant  is  compressed 
and  pumped  out  at  high  speeds  between  the  teeth  and  the  ends.  See  Figure  1.1. 
The  result  of  the  pumping  action  is  noise  and  possibly  thermal  failures.  If  the 
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driver 


Figure  1,1 


Air  and  Lubricant  Flow.  Air  and  lubricant  are  compressed 
and  pumped  out  from  between  the  teeth 
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air  velocity  between  the  teeth  should  reach  sonic  levels,  a  shock  wave  can  form 
which  emajiates  noise.  Drago  [l]  mentions  this  problem  for  high  speed  gears  with 
very  wide  face-widths.  Rosen  [7]  presented  air  flow  noise  as  part  of  an  overall 
examination  of  gear  noise  sources.  In  his  report,  he  computed  the  velocity  of  air 
flow  in  spur  gears  using  incompressible  flow  theory.  By  drawing  curves  of  the 
meshing  volume  and  discharge  area  with  respect  to  mesh  position,  he  generated 
a  graphical  representation  of  the  velocity,  V  =  where  A  is  the  discharge 

area  and  dV  J dt  is  the  time  rate  of  volume  change.  His  analytical  results  indicate 
the  air  velocity  approaches  sonic  levels  for  a  particular  gear  operating  at  20,000 
RPM  and  this  corresponded  to  the  experimental  result  of  a  sudden  rise  in  noise. 
Dudley  recognized  the  thermal  complications  associated  with  fluid  being  expelled 
between  meshing  teeth.  He  attributes  overheating  failures  of  gear  box  and  gear 
teeth  to  high-speed  air  flow  between  teeth.  For  high-speed  gears,  the  gear  teeth 
act  like  blades  on  a  centrifugal  compressor  which  heat  up  due  to  the  friction  with 
the  surrounding  air.  Further  heating  is  caused  by  the  pumping  air  each  time  a 
gear  meshes.  Typical  thermal  failures  are  softening  of  teeth,  scoring  or  pitting 
due  to  oil-film  breakdown  on  overheated  teeth,  local  overloading  of  teeth  due  to 
thermal  distortion  of  tooth-contzxt  pattern,  misalignment  of  gears  and  bearings 
from  distortion  of  the  casing,  and  failure  of  bearings  due  to  overheating.  To 
prevent  problems  of  this  nature,  Dudley  recommends  not  using  spur  and  straight 
bevel  teeth  when  the  pitch  line  velocity  exceeds  10,000  fpm,  because  they  do 
the  poorest  job  of  expelling  air  efficiently.  As  a  guide  to  designers,  he  also  lists 
typical  axial  meshing  velocities  for  helical  gears  and  the  corresponding  severity  of 
thermal  problems  resulting  from  air  flow  [8|.  Buckingham  [9]  notes  that  excessive 
oil  at  the  tooth  mesh  creates  heating  problems  from  the  high  discharge  rate  of 
oil  at  the  mesli.  Smith  mentions  that  noise  can  occur  if  oil  is  trapped  in  the 
roots  of  meshing  spur  gear  teeth  of  high  face  width.  If  oil  cannot  escape  through 


4 


the  ends  of  the  teeth,  high  pressure  develops  and  it  forces  the  gears  apart.  The 
pressurized  oil  shoots  outward  against  the  gear  casing  resulting  in  vibrations  [lO], 
The  purpose  of  this  thesis  was  to  assess  the  magnitude  of  the  fluid  velocity 
that  results  from  meshing  spur  gesj'  teeth  and  to  determine  whether  it  reaches 
levels  capable  of  causing  noise.  The  fluid  velocity  was  analyzed  by  assuming  one- 
dirnensional  flow  and  applying  incompressible  and  compressible  flow  theories. 
It  is  expected  that  this  worl:  will  form  the  basis  for  future  experimental  and 
theoretical  analysis  of  the  fluid  pumping  problem.  Eventually,  this  work  should 
lead  to  a  better  understanding  of  the  hydrodynamics  associated  with  the  meshing 
action  and  permit  the  prediction  of  conditions  causing  noise  and  tooth  failures. 
Predictive  knowledge  can  aid  designers  in  making  more  quantified  decisions  on 
choice  of  geometry  and  pitch  line  velocity. 

A  brief  description  of  the  tooth  geometry,  the  meshing  action,  air  flow  induced 
noise,  and  formation  of  a  shock  wave  is  given.  For  gears  to  run  smoothly  and 
transmit  power  effectively,  the  ratio  of  the  angular  velocity  of  the  driving  gear 
compeured  to  the  driven  must  remain  constant.  The  law  of  gearing  states  that 
the  lines  of  action  at  every  point  of  contact  must  pass  through  the  pitch  point. 
The  involute  profile  is  used  as  a  tooth  proflle  because  its  shape  satisfies  these 
requirements  [11].  The  inv  jlute  profile  is  generated  by  unwinding  a  string  from 
the  circumference  of  a  circle.  See  Figure  1.2.  Definitions  for  the  parts  of  a  gear 
tooth  are  shown  in  Figure  1.3.  The  diametral  pitch  is  the  number  of  teeth  on 
a  gear  per  inch  of  pitch  diameter.  The  pitch  circle  diameter  is  defined  by  the 
number  of  teeth  divided  by  the  diametral  pitch.  The  addendum  is  the  radial 
distance  from  the  pitch  circle  to  the  top  of  the  gears.  The  dedendum  is  the 
radial  distance  from  the  pitch  circle  to  the  bottom  land.  The  base  circle  is  the 
point  where  the  involute  profile  begins. 
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STRING 


Figure  1.2  Generation  of  an  Involute.  The  involute  profile  is  generated 

by  'inwinding  a  taut  string  from  the  circumference  of  a  circle. 
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The  meshing  action  (Figure  1.4)  is  split  into  two  parts;  the  approach  and 
the  recession.  The  aic  of  approach  begins  with  the  initial  contact  of  the  driving 
gear,  Ti  at  the  tooth  tip  of  the  driven,  T2  (solid  lines).  The  contact  slides  up 
to  the  pitch  point  where  the  approach  ends  (dashed  lines).  The  arc  of  recession 
begins  at  the  pitch  point  (solid  lines)  and  ends  when  the  tip  of  the  driving  gear 
loses  contact  with  the  driven  gear  (dashed  lines).  All  of  the  contact  points  fall 
along  the  pressure  line. 

Air  flow  noise  results  from  three  sources;  monopole,  dipole,  and  quadrupole. 
A  monopole  source  occurs  when  air  is  interrupted  at  a  regular  rate.  An  example 
is  a  siren.  A  dipole  source  is  one  in  which  a  moving  stream  strikes  a  solid  object 

and  a  quadrupole  source  is  a  result  of  a  gas  jet.  An  example  of  dipole  noise  is 

_ 

the  aerodynamic  noise  generated  by  fan  blades  moving  thiuugh  air.  Quadrupole 
sources  01  noise  arise  from  a  gas  jet  [12j.  Gas  jet  noise  is  distributed  downstream 
from  the  nozzle  v/ith  high  frequencies  near  the  nozzle  and  lower  frequencies 
downstream.  When  the  critical  pressure  for  a  gas  jet  is  reached,  the  jet  is  choked 
and  a  shock  formation  occurs  which  is  an  additional  source  of  noise  [13].  The 
motion  of  high  speed  gears  probably  gives  rise  to  dipole  and  quadrupole  sources 
of  noise.  The  dipole  source  is  a  result  of  the  gear  teeth  acting  as  blades  moving 
through  air  and  the  quadrupole  source  is  due  to  the  jet  stream  of  air  exhausted 
from  between  meshing  teeth.  Calculations  in  this  thesis  have  shown  that  for 
particular  geometries  and  pitch  line  velocities,  the  jet  velocities  can  reach  sonic 
levels. 

Shock  waves  represent  rapid  changes  in  the  velocity  and  pressure  of  the  flow 
and  tire  modeled  as  a  discontinuity  [14|.  Shocks  are  formed  when  compression 
waves  become  steeper  tind  form  a  discontinuity.  In  the  particular  case  of  meshing 
gears,  when  the  pressure  in  the  mesh  region  increases  enough  for  the  pressure 
ratio,  pressui '  inside  divided  by  ambient,  to  reach  the  critical  ratio  the  jet 
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ARC  OF 
RECESSION 


Figure  1.4  Arc  of  Approach  and  Recession.  The  arc  of  approach  and 
recession  are  shown  for  the  driving  tooth,  Ti,  and 
the  driven  tooth,  Tj.  The  solid  lines  indicate  the  beginning 
of  the  arc  and  the  dashed  lines  the  end  of  the  arc. 
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becomes  choked.  The  flow  expands  downstream  from  the  nozzle  at  supersonic 
velocity  and  a  shock  is  formed  outside  the  nozzle. 

Chapter  2  examines  incompressible  flow  for  the  meshing  spur  gear  teeth  in 
a  stationary  rack  and  gear  using  a  one-dimensional  flow  model.  This  model  is 
most  appropriate  for  heavily  lubricated  gears  in  which  the  composition  of  the 
fluid  trapped  between  the  gear  teeth  is  the  lubrication  oil.  The  assumption  of 
an  ideal  fluid  allows  the  fluid  to  be  considered  frictionless  and  incompressible. 
The  equation  for  incompressible  velocity  is  derived  by  applying  the  continuity 
equation  to  a  control  volume.  The  remainder  of  the  chapter  explains  the 
evaluation  of  the  fluid  discharge  area  and  the  volume. 

A  representation  of  the  flow  problem  for  air  is  presented  in  Chapter  3  which 
accounts  for  the  changes  in  density.  This  model  assumes  the  gears  are  lightly 
lubricated  with  the  primary  discharged  fluid  being  air.  Using  this  assumption, 
the  velocity  was  found  to  be  comparable  in  magnitude  with  the  local  speed  of 
sound.  The  compressible  velocity  equation  was  developed  using  the  First  Law  of 
Thermodynamics  in  addition  to  the  continuity  equation. 

The  hydrodynamic  action  of  the  meshing  gear  teeth  was  simulated  by  a 
computer  algorithm.  These  calculations  gave  information  about  the  following: 
fluid  velocity  at  the  exit  plane  and  the  temperature,  density,  and  pressure  in  the 
tooth  space.  The  results  and  discussion  arc  given  in  Chapter  4.  An  explanation 
of  the  computer  algorithm  and  a  block  diagram  is  presented  in  the  appendix. 


Chapter  2 


INCOMPRESSIBLE  FLOW  THEORY 

2.1.  Introduction 

The  expressions  which  describe  the  time  dependent  velocity  were  developed 
based  on  a  one-dimensional  model  in  which  the  fluid  was  assumed  to  be  friction¬ 
less  and  incompressible.  Thus,  only  the  continuity  equation  and  expressions  for 
the  relevant  geometric  parameters  are  required.  The  model  is  based  on  a  control 
volume  defined  by  the  surfaces  of  the  gear  teeth  and  the  flow  area  at  the  ends  of 
the  gears  and  the  flow  area  defined  by  the  tip  of  the  engaging  tooth  and  the  face 
of  the  corresponding  tooth  on  the  driven  gear.  The  analysis  began  at  the  angle 
of  approach  when  contact  first  occurs.  This  approach  assumed  no  flow  occurred 
before  contact.  This  derivation  is  appropriate  for  a  stationary  rack  and  gear 
and  also  for  a  stationary  ring  and  gear.  However,  for  two  rotating  spur  gears  an 
additional  tangential  velocity  term  for  the  angular  rotation  of  the  control  volume 
should  be  accounted  for  by  the  vector,  rw,  where  r  is  the  pitch  radius  and  W  is 
the  angular  velocity. 

2.2.  Velocity  Equation 

The  velocity  equation  was  derived  by  first  obtaining  expressions  for  the 
control  volume  geometry  and  then  applying  the  continuity  equation  to  solve 
for  the  fluid  velocity. 

2.2.1.  Control  Volume 

As  the  driving  gear  tooth  contacts  a  pinion  tooth,  a  control  volume  is  formed. 
For  the  calculations,  the  driven  gear  was  considered  the  pinion  and  the  driver  the 
gear.  The  control  surface  is  described  by  the  boundaries  of  the  involute  profiles 
of  two  adjacent  pinion  teeth,  their  bottom  land,  and  the  profile  of  the  meshing 
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gear  tooth  and  the  discharge  exit  area.  See  Figure  2.1.  Note  that  the  control 
volume  has  a  variable  geometry.  The  main  effort  in  solving  for  the  velocity  was 
to  obtain  an  expression  for  the  volume  and  flow  area  as  a  function  of  the  relative 
angular  position  of  the  meshing  gear  and  pinion  teeth. 

2.2.2.  Continuity 


The  integral  form  of  the  continuity  equation  states  that  the  time  rate  of 
increzuse  of  mass  in  the  control  volume  is  equal  to  the  net  rate  of  mass  inflow  or 
outflow  to  the  control  volume. 


d 

dt 


ndS, 


(2.1) 


where  p  is  the  density,  V  is  the  volume,  Vr  is  the  velocity  relative  to  the 
exit  plane  of  the  control  volume,  S  is  the  surface  area,  zmd  n  is  a  unit  vector 
normal  to  the  surface.  Assuming  one-dimensional  flow,  an  incompressible  fluid 
of  constant  density,  /),  and  flow  through  a  discharge  area,  A,  gives 


P 


dt 


(2.2) 


Canceling  p  and  evaluating  the  integral  gives 


(2.3) 


Solving  for  the  relative  velocity  gives 


Vr 


A  dt 


(2.4) 


Using  the  chain  rule  gives 


A  dO  dt 


5 


(2.5) 
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CONTROL 

VOLUME 


Figure  2.1  Control  Volume.  A  control  volume  is  formed  as  the  driving  tooth 
contacts  the  driven  tooth. 
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where  dV / d6  is  the  change  in  volume  with  respect  to  position  and  dO J dt  is 
the  time  rate  of  change  in  position.  Defining  the  angular  velocity  of  the  gear, 
(jj  =  dd  fdt,  and  substituting  into  equation  2.5  results  in 


Vr 


1  dV 
A  dQ 


(2.6) 


Dividing  both  sides  by  the  pitch  line  velocity,  UJ  X  rpUcht  gives 

Vr  1  dV 

Vr  =  - =  - rr-  (2-7) 

U)  X  f pitch  A  X  T pitch  do 

This  represents  the  dimensionless  velocity  escaping  through  a  discharge  area,  A. 


2.3.  Geometric  Parzuneters 

The  two  geometric  parameters,  A  and  dV Jd6,  depend  on  the  geometry  of 
the  gears  and  the  position  in  the  mesh  cycle.  The  discharge  area,  A,  was  chosen 
to  have  three  possible  values  depending  on  the  path  of  air  flow:  area  between 
teeth,  end  flow  area,  and  combined  flow  area.  The  rate  of  volume  change,  dV JdO, 
was  calculated  using  a  numerical  derivative  of  equation  2.10.  The  values  for  the 
volume  as  a  function  of  $  were  determined  from  a  2-D  plane  area  which  was 
then  multiplied  by  the  face  width  to  give  a  volume.  Figure  2.2  summarizes  the 
variables  of  the  gear  geometry  used  to  calculate  the  discharge  area  and  volume. 


2.3.1.  Discharge  Area 

Air  escapes  from  the  meshing  region  by  two  paths.  One  path  is  the  direct 
discharge  between  the  meshing  teeth.  The  other  is  flow  out  the  ends  of  the  gear. 
See  Figure  2,3,  The  amount  of  flow  through  either  path  depends  on  the  tooth 
height  and  the  face  width  of  the  gear.  In  the  case  of  wide  gears,  the  discharge 
area  between  the  teeth  dominates  flow.  When  the  face  width  is  narrow,  end  flow 
dominates.  A  combination  flow  results  for  intermediate  face  width  gears. 


END 


Figure  2.3  Flow  Paths  for  Discharged  i  luid.  Air  escapes  either  through 
the  ends  or  between  the  the  teeth* 
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To  determine  the  area  for  discharge  between  the  meshing  teeth,  the  shortest 
chord  distance  between  the  meshing  teeth  is  calculated  and  the  backlash  amount 
is  added  to  the  chord  length.  During  the  approach,  chord2  is  defined  as  the 
distance  between  the  involute  surface  of  the  upper  pinion  tooth  and  the  near  part 
of  the  tip  of  the  engaging  gear  tooth  (trailing  edge  involute,  r,nt;2)'  See  Figure 
2.4.  Thus,  the  area  between  the  teeth  for  flow  is  {backlash  +  chordi)  X 
face  width.  This  is  not  precisely  true  near  the  end  of  the  mesh  cycle  because 
the  involute  portion  of  the  tooth  is  closer  to  the  meshing  tooth  than  the  tip,  but 
at  that  time  in  the  mesh  the  chord  is  closed  off  and  the  distance  is  reasonably 
represented  by  the  backlash  amount.  The  backlash  is  the  amount  the  teeth  are 
cut  thinner  than  the  theoretical  dimensions.  The  backlash  accommodates  any 
expansion  of  the  teeth  due  to  temperature  effects. 

The  end  flow  area  for  one  end  equals  the  control  volume  divided  by  the  face 
width.  The  total  end  flow  area  is  twice  the  single  end  area  and  is  given  by 

total  end  area  =  2  x  {volume/  face  width).  (2.8) 

The  combined  area  was  chosen  to  be  the  sum  of  the  area  between  the  meshing 
teeth  and  the  end  flow  with  a  factor  related  to  the  tooth  height  and  face  width. 
The  combined  area  is  given  by 

combined  area  =  between  teeth  +  (end  area  X  /)  \  ^ 

'  ^  ^  face  width  ' 

(2.9) 

where  f  is  the  percentage  of  end  area  used.  This  equation  wzis  chosen  as  a 
possible  representation  of  the  distribution.  It  is  based  on  the  idea  that  the  fluid 
will  tend  to  seek  the  shortest  path  of  least  resistance.  Using  the  ratio  of  the 
tooth  height  to  the  face  width  resembles  the  length  of  path  the  fluid  may  travel. 
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Figure  2.4  Chord  Between  the  Teeth.  The  shortest  chord  length  between 
the  teeth  was  used  to  calculate  the  discharge  area 
between  the  teeth. 
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2.3.2.  Volume 

The  rate  of  volume  change  with  respect  to  position,  9,  was  calculated  by 
taking  a  numerical  derivative  of  equation  2.10.  The  volume  was  described 
previously  as  being  bounded  by  two  pinion  teeth  and  the  meshing  gear  tooth. 
The  volume  was  calculated  by  multiplying  the  two-dimensional  plane  area,  Figure 
2.5,  by  the  constant  face  width  of  the  gear. 

volume  =  {2  — D  planearea)  X  /ace  width.  (2.10) 

The  2  —  D  planearea  is  the  sum  of  the  area  between  two  pinion  teeth,  Bareat 
minus  the  area  overlapped  by  the  meshing  gear  tooth,  overlaparsa^  and  minus 
the  section  no  longer  part  of  the  control  volume,  excludedarea-  A-  description 
follows  for  each  component  of  the  2  ~  D  area. 

2  D  planedYgd  —  Ba^^d  ~  overlap (itrol  excludeduTcci’ 


The  area  between  two  pinion  teeth  is 

D  ^  (fap)^  -  TT  (rdp)2  -  N  X  tooth 

N 


area 


'area 


(2.12) 


where  N  is  the  number  of  teeth  and  the  area  of  a  pinion  tooth  is  given  by 
equation  2.13.  Figure  2.6  gives  a  description  of  the  limits  for  the  tooth  area 
calculation. 


toOthdTca, 


where  tr^^g 
pinion,  rinv 


tr^p  Tap 

j  j  r  dr  d(f) 

0  Tbp 


+ 


tT*6p 

/  r  dr  d(f) 

0  r^p 


Ep  Tap 


(2.13) 


0  Tinv 

is  the  thickness  in  radians  of  the  tooth  at  the  base  circle  of  the 
is  an  approximation  for  an  involute  radius,  and  Ep  represents  the 
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■  OVERLAP  ARE-^ 


Figure  2.5  2-D  Plane  area.  The  volume  was  determined  by  multiplying  the 

2-D  plane  area  by  the  gear  face  width. 
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arc  in  radians  from  the  beginning  of  the  tooth  profile  at  the  base  circle  to  the 
tooth  tip  [9]. 

Vbp  ^bp 

Evaluating  the  integrals  for  tootharea  gives 

tootharea  —  ~  ^bp)^^bp  +  2^^ ~  ^dp)^^bp 

-  [rlpEp  -  rlp{\  xEp  +  3^/®  X  (2.15) 

0 

+  i  X  X  El/% 

The  area  overlapped  by  the  meshing  gear  tooth,  overlapareaz  is  approx¬ 
imated  by  integrating  the  area  of  the  gear  tooth  profile  inside  the  addendum 
of  the  pinion.  The  integrals  used  depended  on  the  location  of  the  driving  gear 
tooth  centerline  with  respect  to  the  x-axis  and  whether  or  not  the  base  circle 
radius  of  the  gear  overlapped  the  addendum  circle  of  the  pinion.  See  Figure  2.7. 
If  the  base  circle  radius  overlapped  the  pinion  addendum  during  the  mesh,  four 
sets  of  integrals  resulted:  Sets  1  through  4.  In  the  case  of  the  bzise  circle  of  the 
gear  not  ovcrla^y—  ^  the  addendum  circle  of  the  pinion,  only  two  sets  of  integrals 
resulted:  Sets  5  and  6.  For  the  calculation  of  over/apcreo?  approximations  for 
the  involute  profile  of  the  gear  and  the  equation  of  the  pinion  addendum  circle 
were  developed. 

The  involute  profiles  for  the  driving  gear  tooth  are  the  leading  side  finvi, 
and  the  trailing  side,  rinv2  see  Figure  2.8.  The  profiles  were  approximated  as 

r.nvi  =  rbgil  +  (2.16) 

and 

^inv2  “  ^bg(^  ^))  ^  • 


(2,17) 


PINION 

ADDENDUM 

CIRCLE 


Base  circle  overlaps  addendum.  When  the  base  circle  of  the  gear 
falls  inside  the  addendum  of  the  pinion,  four  sets  of  area  integrals  are 
required. 


Base  circle  outside  of  addendum.  When  the  base  circle  of  the  gear 
remains  outside  of  the  pinion  addendum,  two  sets  of  area  integrals  are 
required. 


Figure  2.7  Determining  Integral  Sets. 
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The  approximations  were  developed  using  Buckingham’s  [9]  vector  representation 
of  <j>.  See  Figure  2.8. 


4 


—  r* 

'  invl  bg 


Tbg 


—  arctan 


63 


(2.18) 


Defining  gives 


<j>  =  \/f^  —  1  —  arctan  \/f^  —  1. 


(2.19) 


Expressing  arctan  —  1  as  an  infinite  series  and  eliminating  higher  order 
terms  gives 


(2.20) 


and  solving  for  r^^vi  gives 


Unvi  =  rbg{l  +  (30)"/^)'/".  (2.21) 


A  binomial  expansion  and  elimination  of  higher  order  terms  allows  even  a  more 
convenient  form  for  integrations. 

rinvi  +  (2.22) 

A  comparison  of  the  actual  involute  to  the  approximation  for  one  set  of  gears  is 
given  in  Figure  2.9.  The  other  limit  of  integration, rp,>i,  is  the  addendum  of  the 
pinion  with  respect  to  the  origin  of  the  gear.  It  is  derived  in  Figure  2.10  by  using 
a  coordinate  transformation. 

^pin  =  pg  "h  ^pp)  9  VQ  (2.23) 

The  positions  of  the  gear  tooth  centerline  where  9  is  defined  from  the  X  — 
axis  for  the  beginning  and  end  points  for  each  set  of  integrals  were  designated 
respectively  as  stavti  and  stopi  for  Set  1  and  start2  and  Gtop2  for  Set  2  etc. 
When  the  angle  of  the  centerline  of  the  gear  tooth  falls  between  angle  starti 


Figure  2.9  Comparison  of  the  actual  radius  profile  to  the  approximate  profile. 
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Figure  2.10  Addendum  of  the  Pinion.  The  addendum  of  the  pinion, 
Tpin,  using  a  coordinate  transformation  with  respect  to 
the  center  of  the  driving  gear,  (0,0). 
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and  stopi  integral  Set  1  is  used.  The  overlap  was  calculated  by  first  integrating 
the  complete  area  under  the  tooth  profile  and  then  subtracting  the  area  under 
the  pinion  addendum  curve.  The  tooth  profile  uses  a  local  coordinate  system  for 
the  angle  integration  d(i>  and  the  pinion  addendum  uses  a  global  system  dO.  See 
Figure  2.11  for  details. 

The  beginning  of  Set  1  occurs  when  the  tooth  tip  of  the  leading  edge  contacts 
the  addendum  circle  of  pinion. 


start  I  =  arccos 


' ap  '  ag 


(rpg  -t-  Tpp) 


2-, 


+  (2.24) 

^  Tog 


2rag{rpg  +  Tpp) 

where  tpy  is  the  thickness  of  the  gear  at  the  addendum  circle.  Stopy  ends  the 
set  when  the  trailing  edge  involute  enters  addendum  circle.  Figure  2.12  gives 
details  of  angles. 

tna  , 

(2.25J 


stopi  =  starti - 

Tag 


The  integrals  for  Set  1  are 

C~X  Tintil 


C~X  rintil  D  —  X  1‘ag  D  ^pin 

overlaparea  =  j  j  r  dr  d(l>  +  j  j  r  dr  d(j)  —  j  j  rdrd$, 


B-X 


C-X  0 


(2.26) 

where  angle  X.  is  from  the  Xaxia  to  leading  edge  of  the  gear,  angle  B  is  the 
point  on  the  involute  profile  that  intersects  the  addendum  of  the  pinion,  angle 
C  is  the  end  of  the  involute  profile  and  angle  D  is  the  point  of  intersection  for 
the  two  addendum  circles. 


Evaluating  each  integral  over  its  limits  gives: 


j  j  r  dr  =  ^rl,[{C  -  X)  -  {B  -  X) 

B-X  0 

+  -  X  3^/^  X  ((C  -  -{B-  X)®/'') 

5 

+  ^  X  3^/"  X  ((C  -  X)"/"  -{B-  X)"/')] 


(2.27) 
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The  beginning  of  Set  1  occurs  when  the 
tooth  tip  of  the  leading  edge  contacts  the 
addendum  of  the  pinion 


Set  1  ends  when  the  trailing  edge 
involute  enters  the  addendum 


circle,  r^_ 

dP 


SAMPLE  INTEGRAL  FOR  SET  1 


Figure  2.12  Integral  Set  1. 
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and 

D  —  X  ^ag 


j  j  rdrd4  =  \rl^[(D-X)-(C-X)\ 


C-X  0 

and 
D  ^pin 


j  J  r  dr  dd  =  ~(sin  2D  —  sin  2B){rpg  +  Tpp)^  +  f‘ap{D 

2(rpg  +  rpp')T0_p 


-B) 


B  0 


L2  V  rap  ' 


- -sinBjl  - 

2  ^  rap 

1  rap  /  .  ffrpg  +  rpp 

- (arcsinl  ( - )  smU  ) 

2  Tpj  +  Tpp  Tap 

—  arcsmf(-^^ - ^)sin5^)  . 

^oo  '  - 


(2.28) 


(2.29) 


The  beginning  of  Set  2  is  at  the  end  of  Set  1. 


start2  =  stopi  . 


(2.30) 


Stop2  occurs  when  the  base  circle  at  the  leading  edge  of  involute  enters  the 
addendum  of  pinion.  See  Figure  2.13. 

rip  ~  rig  —  {Tpg  +  Tpp)' 

-2r65,(rpg +  rpp) 


stop2  =  arccos 


2rr 


ag 


The  integrals  for  Set  2 

C  —  X  Ttinvi 


O — j\.  Ttlnvi  ^  ^  ^ag 

erlaparea  =  +  j  j  r  dr  dcj)  +  I  j  r  dr  d(f) 


ov 


B-X  0 


C-X  0 
G  ^ pit 


I  rdrd(l)  —  j  j  r  dr  d9, 


F-X  0 


B  0 


(2.31) 


(2.32) 
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where  angles  X,  C  are  as  defined  for  Set  1  and  angle  F  is  from  the 
X  —  axis  to  the  start  of  the  trailing  edge  involute  and  angle  G  is  the  intersection 
point  for  the  trailing  edge  involute  and  the  addendum  circle  of  the  pinion. 

Evaluating  each  integral  over  its  limits  gives 


j  j  r<ird^  =  ir’,[(C-X)-(B-X) 

B-X  0 

+  -  X  3"/^  X  ((C  -  -{B-  X)5/=*) 

5 

+  —  X  3"/"  X  {{C  -  -{B-  X)"/")] 

28  ^ 

and 

F-X  rag 

I  J  r  dr  d<l.  =  ~rl,  [(F  -  X)  -  (C  -  X)]  (2.34) 

C-X  0 

and 


J  j  r  dr  d^  =  ^-rl,  [(tr,,,  -  (G  -  X))  -  -(F-X)) 

F-X  0 

+  J  X  3'/“  X  ((trt,  -  (G  -  X))=/” 

-  (tn,  -(F-  X))‘'^) 

+  -  X  3^/’  X  {{tn,  -  (G  -  X))"/” 

28 

-  (tri,,  -(F-  X))'/’)] 


(2.35) 
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and 


G  •'pin 


j  J  r  dr  dd  = -{sin  2G  —  sin  2 B){rpg Tpp)^  +  r^p{G  —  B) 

2(rpj  +  r pp)r ap 


B  0 


-  sin  Gi  1  —  ( 

L2  V  r 


ap 


- -sinBjl-i^^^'^^^^ysin^B 
2  \  r^p 


ap 


2rpg  +  rpp 


(arcsinf(— ^ - —)  sin 

rap 


.  f ,^pg  +  .  d'\'4 

—  arcsmn - )sin±>jj  . 

^ap 


(2.36) 


Set  3  starts  with  the  end  of  Set  2. 


starts  =  stops- 


(2.37) 


Stops  occurs  when  the  base  circle  at  the  trailing  edge  of  involute  enters  the 
addendum  of  pinion.  See  Figure  2.14. 


stops  -  -Harts  —  trtg-  (2.38) 

where  trig  is  the  thickness  of  the  gear  at  the  base  circle  in  radians.  Set  3  has  two 
cases  for  its  integral  sets.  Case  one  the  leading  edge  of  the  gear  tooth  is  above 
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Set  3  ends  when  the  base  I 
circle  at  the  trailing  edge  involute 
enters  the  addendum  circle,  r^p 


The  edge  of  the  tooth  is  above 


The  edge  of  the  tooth  is  below 
the  x-3xis 


Figure  2.14  Integral  Set  3. 
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the  x-axis  and  Case  two  it  is  below  the  x-axis.  The  integrals  for  calculation  are 

C  —  X  Ttnvl  F  —  X  *"fag 

overlaparea  =  j  J  T  dr  +  j  j  r  dr  d<l> 

c-x  0 


0  0 
G  — ^  r^nva 


+  j  j  r  dr  d(f) 


F  -X  0 


case  one 


G  ^pin 

—  J  j  r  dr  do 

X  0 


(2.39) 


case  two 


G  ^pin  0  ^pin 

-  j  j  r  dr  dO  -  j  J  r  dr  dO, 

0  0  X  0 


where  angles  X,  C,  G  have  been  defined  previously.  Evaluating  each  of 
the  integrals  gives 


C  —  X  r^nvl 

j  j  rdrd^  =  ir2,[(C-X)  +  jX3‘/^x  ((C-X)'/=) 

0  0 

+  ^  X  3*^’  X  ((C  -  X)'/")] 

F-X 

j  j  r<ir<f^  =  iri,[(f-X)-(C-X)] 
c-x  0 


(2.40) 


(2,41) 
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G — X  r^nvs 


F-X  0 


I  J  rdrdij)  =  [(trjj  -  (G  -  Jf ))  -  (irt,  —  {F  —  X)) 

+  3^/’  X  5  X  ((trt,  -  (G  - 

-  (Jrt,  -  (X  -  X))'/“) 

+  -  X  3^“  X  {(tn,  -  (G  -  X))'/^- 
28 

(tri,  -  (F  -  X))’/")] 


G  >'ptn 


X  0 


J  J  r  dr  dO  = -{sm2G  —  sm2X){rpg  +  rpp)^  +  r^p{G  —  X) 

2(rpg  +  Tpp)^ap 

i sin G, /i ■•^-£E)^sin‘G 

2  V  Top 


- -sinX./l '^’'-'’g)^sin^X 

2  V  '■«P 

1  ^PP\ 

- (arcsml  ( - j  sinLr  1 

r«o  ^  r  ap 


2  Tpg  +  /  pp 

Tort  +  Tt 


—  arcsinf(-^^ - sinX^) 

^ap 


G  r 


0  0 


j  J  r  dr  d9  =  ^(sm2G)(rpg  +  rpp)^  +  ^op(^) 

—  2(rpg  +  ^pp)^ap _ _ 

i  sin  G. 


ai7 


1  r 


+  -■ 


ap 


2  r  pg  +  r  pp 


(arcsin  f  ( - —)  sin  ) 

'  Tap 


(2.42) 


(2.43) 


(2.44) 
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0 


j  j  r  dr  d9  =  -^(sin2X)(rpg  +  Tpp)^  -  rJp(X) 

+  2(rpj  +  rpp)rop 


X  0 


rpjjtlPPy  3i^2  X 


(2.45) 


ap 


1  r 


ap 


H - 

2  'f* pg  “H  ^ pp 

Set  4  staxts  at  the  end  of  Set  3. 


f  ap 


starts  =  stop^.  (2.46) 

Set  4  stops  when  the  centerline  of  the  gear  tooth  reaches  the  x-axis.  See  Figure 
2.15. 

Stop^  —  0.0.  (2-47) 

The  integrals  are 


C  —  X  r^nvi 


F  —  X  ^ag 


overlaparea  =  J  j  r  dr  d(f)  +  j  j  r  dr  d<f) 


0  0 
H  ^  ^inv3 


C-X  0 
H 


+  j  j  r  dr  d(l)  ~  j  j  r  dr  dO 


(2.48) 


F-X  0 

0  ^p*n 


0  0 


—  J  f  r  dr  d9, 


X  0 

where  angles  X,  C,  &c  F  have  been  defined  previously  and  angle  H  is  the  angles 
from  the  X  —  axis  to  the  trailing  edge  of  the  driving  tooth.  Evaluating  each 
integral  over  its  limits  gives 

3 


C  X 


I  I  r  dr  d<p  =  [{C  -X)  +  3^/^  X  -  X  ((C  -  X)^/^) 

+  -  X  3^/"  X  ((C  -  X)"/")] 


(2.49) 


The  beginninij  of  Set  4  occurs 
at  the  end  of  Set  3 


Figure  2.15  Integral  Set  4 
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F-X  r, 


r  — ^  lag 

I  I  rdrd<j>^^-rl,[{F-X)-{C-X)] 


(2.50) 


C-X  0 
H  X  Tjnva 


j  J  r  dr  d(p  =  [(tri^  -  {H  -  X))  -  [tr^g  -  {F  - 


X)) 


F-X  0 


+  X  I  X  {(trig  -(H  -  X)fl^ 

-  (tn,  -(F-  X)fi^) 

+  ^  X  3^^  X  {(trig  -  (ff  -  X)yl^ 

-  (trig  -(F-  X))'/=)] 


0  0 


j  J  r  dr  d6  —  ^^{sm2H){rpg  +  rpp)  +  r^p(JJ) 

~  2(rpg  +  Tpp) 

_ 

2  V  rap 


+  - 


ap 


—  (ajcsin ( ( - —)  sin H^) 

PP  ^ap 


2rpgB  +  r 


(2.51) 


(2.52) 


0  ^ptrv 


j  j  rdrd9  =  -^(sin2X)(rpg  +  fpp)"  -  rlp{X) 


X  0 


+  2{rpg  +  rpp)rap 


-sinxjl  -  (^Pg'^^PP)2sin2X 

L2  V  r, 


ap 


,  ^  P“p  /  •  (t^pg'^^pp\  .  v'^^ 

9r  (arcsm(^( - - - )  smXj) 

^»pg"r'pp  'ap  J 


(2.53) 


When  the  bcise  circle  of  the  gear  does  not  overlap  the  pinion  addendum  only 
two  sets  of  integrals  result:  Sets  5  and  6.  Set  5  begins  when  the  leading  edge  of 
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the  tooth  tip  enters  the  addendum  circle  of  the  pinion.  See  Figure  2.16. 


etart^  =  arccos 


op 


^ag  i^pg  +  ^pp) 


2-, 


2r 


(2.54) 


ag 


-2rag{rpg  +  rpp.) 

Set  5  ends  when  the  trailing  edge  involute  enters  the  addendum  of  the  pinion 


stops  =  starts - 


T*  It 


^ag 
^ag 

U —  R  ^ ag 


(2.55) 


overlaparea  =  j  J  T  dr  d(t>  +  j  j  r  dr  d(l) 


S-R  0 

U  ^pin 


T-R  0 


(2.56) 


—  j  r  (ir  dd, 


where  angle  R  is  from  the  x-axis  to  the  leading  edge  of  the  driving  tooth,  angle 
S  to  the  intersection  of  the  leading  involute  profile  and  the  addendum  of  the 
pinion,  angle  T  to  the  end  of  leading  profile,  angle  V  to  the  beginning  of  the 
trailing  edge  profile,  and  angle  W  to  the  intersection  of  the  trailing  edge  profile 
and  the  addendum  circle  of  the  pinion.  Evaluating  each  integral  gives 
T  It  r  *  n  i>  '■ 

J  j  r dr d4  =  ~rl[(T  -  R)  -  l,S  -  R) 


S~R 


+  3=/"  X  -  X  ((T  -  -  (5  - 

5 

+  -  X  X  ((T  -  Rf^^  -  (5  -  RY^^)] 
28  ^ 


(2.57) 


U^R 


f  J  rdrd,i,^.^-rl,\{U-R)-(r-R)] 


T-R  0 


(2.58) 
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U  ^*p%n 


j  J  r  dr  d9  =  ^{sm2U  —  sin  2 S){rpg  +  rpp)^  +  r‘^p{U 

—  2(rpg  +  rpp)rap 


-S) 


s  0 


-smuJl- 

2  V 


-ls\nSJl-C~-^)^sin'^S 


ap 


1  r. 


ap  [  .  f ,^pg  +  ^pp 

- larcsmi  ( - 

2  r pg  "h  r pij  ''  p 


'  pg  ^  '  P'3 

Jvc  +  r.. 


ap 


)  sin 


—  arcsin((-^^' - — )  sin /S')) 

^  rap  M 


The  beginning  of  Set  6  starts  at  the  end  of  Set  5. 


starts  =  stops. 


(2.59) 


(2.60) 


Set  6  ends  when  the  centerline  of  the  gear  coincides  with  the  x-ajcis.  See  Figure 
2.17. 


stops  =  0-0  • 

T-R  Vinvi  V-R  rag 

overlaparca  =  j  j  r  dr  d<l>  +  j  j  r  dr  d(f) 


S~R  0 

W^  —  R  Tinvl 


T-R  0 


4-  j  j  r  dr  d<l) 


for  case  one 


V-R  0 

W  rpin 


for  case  two 


—  j  J  r  dr  d9 

R  0 

W  ^ptn  0  ^pin 

—  j  j  r  dr  d9  —  j  j  r  dr  dO. 


0  0 


R  0 


(2.61) 


(2.62) 
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Evaluating  each  integral  gives 


T—R  finvl 

j  J  rdrd(f>  =  -rig  [(T  -  R)  -  {S  -  R) 

S-R  0 

+  -  X  3"/"  X  ((T  -  -  (S  - 

5 

+  ^  X  3"/"  X  ((T  -  i?)"/"  -{S  -  Ry^^)] 

V-i^  rag 

I  j  r  dr  =  ^rl^[(V  -  R)  -  [T  -  R)] 

T-R  0 


(2.63) 


(2.64) 


W  — 

j  j  rdrd4>  =  [{trbg  -  {W  -  R)]  -  [tr^g  -  (V  -  R)) 

V-R  0 

H-  ^  X  3=/=*  X  {{trtg  -{W-  R)y^^ 

-  {trtg  -{V-  R)y^^) 

+  -  X  3^/=*  X  {{trbg  -{W~  R)y^^ 

28 

-  (irt,  -  (V  - 

(2,65) 

W  ^pirv 

j  j  r  dr  do  = -{sin  2W  —  sm2R){r pg  +  Tpp)^  +  rlp{W  —  R) 

R  0 


^{^pg  +  1"pp)'f‘ap 


-sinvyji  -  {-^^^^ysWW 

L2  V  ^ap 


--sinig,/l-(^P^  ^-gP-)^sin^i^ 

2  V  ’’“P 

_  1  ’•-P  (.rcsin(('''’^'*'-^^-)sin>y) 

2  Ppg  4"  r pp  r fxp 

.  +  ^PP^  D^^ 

—  arcsinl  ( - )  sm  it  j  j 

an 


(2.66) 
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W  rpin 


I  j  r  dr  d9=  ^(sin  2W){rpg  +  rpp)^  +  rlp{W) 

0  0 

-  2{rpg  +  rpp)rap 


V9 


+  r 


ap 


singly 

pg  +  Tpp^ 


(2.67) 


_l - - (arcsin  ( ( - —)  sin  ) 

^  ^pg  “I"  ^PP  ^ dp 


0  ^ptn 


j  J  r  dr  d9  =  —-{sin2R){rpg  +  rpp)^  —  rlp{R) 
R  0 

+  2{rpg  +  r  pp)r  dp 


-sini2\/l-( 


!Vi±rffi)2  R 


(2.68) 


ap 


+ 


op 


2r, 


■  (arcsin — —)  sin  R^) 

^ap 


pg  -r  fpp 

The  area  subtracted  from  the  space  between  two  pinion  teeth  is  derived  from 
the  dischairge  area  and  is  designated  as  the  excludedarea-  When  the  chord 
length  along  the  disharge  area  is  known,  the  area  from  the  exit  plane  out  to  the 
addendum  of  the  pinion  can  be  calculated.  The  area  above  the  pinion  involute 
prohle  can  also  be  calculated.  See  Figure  2.18.  At  the  point  in  the  meshing 
cycle  when  the  tooth  tip  of  the  driving  gear  is  no  longer  the  closest  point  to  the 
driven  tooth  (see  chord  2  description),  another  approximation  is  used  to  compute 
the  excluded  area.  A  cubic  equation  is  fit  using  the  last  area  calculated  and  an 
estimation  of  the  final  value.  The  integrals  are 

Pa  ^ ap  ^ap 

excludedarea  =  J  j  r  dr  d(f)  +  II  r  dr  d(j) 


0  Pl 


Pr 


'p 


—  j  j  r  dr  d(f>. 
Pr  0 


(2.69) 


P5  P4  P6  PC  12 
TO  DASHED  LINE 


Figure  2.18  Excluded  Area.  Calculation  of  the  area  excluded  from 
the  coiitrol  volume. 
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Evaluating  ^ach  integral  gives 

Pq  rap 


J  j  r  dr  d(j>  =  "  ^l) 


0  Pl 

Ep  Tap 


j  J  r  dr  d<l)  =  -»‘op(-^p  “ 


Pr  0 


Pr  0 


+  -i  X  3^/^  X  {e;/^  -  p'A]- 
28 


The  total  plane  area  subtracted  equals 

exludedarea  =  “-fsC^ttp  ^l) 


-  Pr) 


i((Pp  -  Pr) 

+  i  X  3=/='  X  (£'/"  -  P"’^) 

+ i  X  3’/=  X  (Ey^  -  P7^“)i. 
28 


(2.70) 


(2,71) 


(2.72) 


(2.73) 


This  chapter  has  presented  a  method  of  solution  of  incompressihle  flow  of 
air  between  spur  gear  teeth.  Through  the  use  of  assumptions,  a  numerical 
representation  of  the  change  in  discharge  area  and  change  in  volume  were  given. 


Chapter  3 


COMPRESSIBLE  FLOW  THEORY 

3.1.  Introduction 

By  taking  into  account  density  changes,  an  approximation  of  the  velocity 
for  a  perfect  gas  resulting  from  the  meshing  of  a  gear  and  stationary  rack  was 
developed.  To  develop  the  velocity  equation,  the  First  Law  of  Thermodynamics 
was  applied  to  relate  the  fluid  properties  along  a  streamline  from  inside  the 
control  volume,  where  the  velocity  is  assumed  to  be  low,  to  the  exit  plane.  To 
determine  the  pressure  relationship,  the  continuity  equation  was  applied  to  the 
control  volume. 

An  ideal  gas  and  one-dimensional  flow  model  with  constant  entropy  and 
adiabatic  flow  was  assumed.  This  gives  constant  properties  for  the  velocity, 
pressure,  density,  and  temperature  across  the  flow  area.  The  assumption  of 
constant  entropy,  ds  =  0,  allows  any  thermodynamic  variable  to  be  found  in 
terms  of  any  other  thermodynamic  variable.  Isentropic  flow  can  be  assumed 
because  friction  effects  are  small  over  the  short  distance  of  accelerating  flow  and 
heat  transfer  is  small  because  of  the  short  time  intervals.  The  analysis  began  at 
the  angle  of  approach  when  contact  first  occurs  between  the  driving  and  driven 
teeth.  Before  contact  occurred,  it  was  assumed  that  the  pressure  inside  the 
control  volume  was  the  same  as  the  ambient  and  no  flow  occurred.  The  discharge 
area  evaluation  and  the  volume  are  the  same  as  the  incompressible  case  and  were 
described  in  Chapter  2. 

3.2.  Velocity  Equation 

Applying  the  First  Law  of  Thermodynamics  along  a  streamline  gives  the 
velocity  equation.  The  pressure  relationship  of  the  velocity  equation  is  derived 
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by  applying  continuity  to  the  control  volume. 

3.2.1.  Energy  Equation 


Applying  the  First  Law  of  Thermodynamics  along  a  streamline  for  adiabatic 
conditions  gives 

L  L 

hi  +  -  =  hg  +  ,  (3*1) 

2^c  ^Qc 

where  Vi  is  the  velocity  inside  the  control  volume,  Vg  is  the  velocity  at  the 
exit  plajie,  hi  and  hg  are  the  enthalpy  inside  the  control  volume  and  at  the 
exit  plane,  gg  is  a  dimensional  constant,  and  Cp  is  the  specific  heat  at  constant 
pressure.  Rewriting  the  first  law  with  =  0  and  h  =  CpT  gives 


CpTi  =  CpTg  H - 

‘^9c 

Solving  for  the  velocity  squared  at  tne  exit  plane  gives 

=  2gcCp{Ti  -  Tg) 


Vg  =  2ggCpTi[l  -  — ) 


(3.2) 

(3.3) 

(3.4) 


Dividing  both  sides  of  the  equation  by  the  pitch  line  velocity,  W  X  Tpg,  and 
defining  Vg  =Vgl[ijJ  X  Tpg)  and  T  =  TifT^Q  results  in 


-2  ^Qc^pl'oo 

vt  = - - - r 

«  /.  .  -  \2 


T 

T(i  -  — i 

(wxr„)=  '  rr 


(3.5) 


Substituting  for  TgfTi  with  the  isentropic  relationship  P  *  where  P  = 
Pi  J Pg  and  taking  the  square  root  of  both  sides  gives 


Vg  = 


oo^. 


w  X  r 


-T(i-P— 


(3.6) 


vg 


The  pressure  variable,  P,  is  evaluated  in  section  3.2.2.  by  applying  the  continuity 
equation  to  the  control  volume.  The  temperature,  T,  is  evaluated  by  using  the 
isentropic  relationship 


Ti  Pi 


.(3-7) 
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The  other  relationship  can  be  written  as 


-  = 

Pt  V. 


(3.8) 


3.2.2.  Pressure  Equation 


(3.9) 


Writing  the  continuity  equation  for  the  control  volume 

d 

I  PidV  + 

’s 

where  Vr  is  the  relative  velocity  to  the  exit  plane,  p,-  is  the  density,  V  is  the 
volume,  and  S  is  the  surface  of  the  control  volume  .  For  he  1-D  case 


^  Jv  Is 


+  PeAVr  =  0. 


Expanding  the  equation  gives 


dpi  dV 


Multiplying  by  gives 

^  dpi  1  dV  pg  A 

Jilt  vlt 

Using  the  isentropic  case  and  making  substitutions  gives 

1  1  dP  I  dV  fPgx\/kA 


+ 


+ 


\p.) 


kP,  dt  V  dt  '  ^Pi 

Defining  P  =  Pi  f  Pg  and  substituting  gives 

\  I  dp  \  dV  ^^i/k  ^ 

Using  the  chain  rule  for  dP fdt  and  dV / dt  gives 

1  1  dP  d9  1  dV  dO  ~  ./u  A 

- L. - - - L  P  - t;  Q 

kp  dO  dt  V  de  dt  y  ^ 


(3.10) 


(3.11) 


(3.12) 


(3.13) 


(3.14) 


(3.15) 
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Solving  for  dPfdO  gives 


dP  r  PdV 
d9  "IV  dd 


p(k-i)/k 

V  u  \  ' 


(3.16) 


The  pressure,  P  in  equation  3.16,  was  solved  using  the  Runge-Kutta  fourth  order 
formula.  The  step  size,  A^,  was  chosen  by  running  the  program  with  a  small 
step  size  and  then  cutting  the  value  in  half  and  running  the  program  again  to 
compare  it  to  the  previous  trial  step  size.  When  the  difference  in  the  results  was 
small,  less  than  \%,  the  step  size  was  set. 

Until  critical  conditions  are  reached,  the  conditions  at  the  exit  plane  must 
be  equal  to  the  ambient  pressure  =  p^Q  and  P^  =  Pfx>-  The  meshing 
teeth  cm  be  seen  to  form  a  converging  nozzle.  See  Figure  3.1.  If  P^  were  larger 
than  Poo,  the  stream  would  expand  laterally  upon  leaving  the  nozzle;  however, 
an  increase  in  area  at  subsonic  speeds  would  cause  the  pressure  to  rise  further. 
Since  Pqo  is  the  ultimate  pressure  reached,  P^  cannot  be  larger  than  foo  (14]. 
Once  the  pressure  in  the  cavity  divided  by  the  ambient  has  reax:hed  the  critical 
ratio,  the  flow  becomes  choked  and  remains  at  Mach  1. 

The  compressible  velocity,  equation  3.6,  approximates  the  condition  of  lightly 
lubricated  gears  where  it  can  be  assumed  that  the  fluid  is  only  air.  By  using  this 
equation,  it  can  be  determined  whether  the  velocity  reaches  levels  capable  of 
producing  shocks  and  causing  noise.  The  true  velocity  should  include  the  effects 
of  the  exit  plane  moving  inward  and  be  accounted  for  in  the  energy  equation. 
By  neglecting  it,  a  difference  on  the  order  of  5%  is  introduced.  Also,  for  two 
rotating  gears  the  effect  of  the  rotating  control  volume  should  be  included  in  the 
energy  equation. 


CONVERGING 

NOZZLE 


Figure  3.1  Converging  Nozzle.  The  meshing  gear  teeth  forming  a  converging 
nozzle. 


Chapter  4 


RESULTS  AND  DISCUSSION 

4.1.  Introduction 

When  the  driving  gear  contacts  the  driven  pinion  tooth,  fluid  begins  to  flow 
out  from  the  volume  between  the  meshing  teeth.  An  estimation  of  the  fluid 
velocity  relative  to  the  exit  plane  was  obtained  by  using  incompressible  and 
compressible  flow  theories.  Figures  4.1  thru  4.7  contain  incompressible  flow 
results  and  Figures  4.8  thru  4.24  have  the  compressible  flow  results.  Figures 
4.25  through  4.27  contain  comparisons  of  incompressible  and  compressible  flow 
for  each  of  the  assumed  discharge  areas.  The  geometric  parameters  for  the 
incompressible  flow  figures  are  given  in  Table  1.  Table  2  contains  the  compressible 
geometry  parameters  and  Table  3  has  the  comparison  geometries. 

It  was  assumed  no  flow  occurred  before  the  driving  tooth  contacted  the 
driven  tooth  and  no  thermal  expansion  of  the  teeth  due  to  temperature  rises. 
The  percent  of  mesh  axis  on  the  plots  indicates  the  position  of  the  driving 
tooth  as  a  percentage  of  the  arc  from  the  initial  contact,  zero  percent  mesh, 
to  full  mesh  position,  100  percent  mesh.  The  velocity  is  given  as  a  ratio  of 
the  computed  fluid  velocity  to  the  pitcl  line  velocity.  The  three  discharge  areas 
for  analysis  were  end  flow  area,  teeth  flow  area,  and  combined  flow  area.  The 
endflow  Wcis  chosen  to  represent  the  meshing  of  narrow  face  width  gears.  For 
this  case,  the  fluid  is  expected  to  flow  primarily  out  the  end  openings  of  teeth. 
The  other  extreme  is  wide  face  gears.  A  wide  face  width  should  tend  to  force 
the  discharged  air  out  between  the  teeth  and  Wcis  designated  teeth  flow.  The 
combined  flow  uses  equation  2.9  to  account  for  a  combination  of  the  two  flow  areas 
by  using  the  height  of  a  tooth  and  the  face  width.  The  influence  of  the  geometric 
parameters  including  backlash,  diametral  pitch,  drive  ratio,  and  pressure  angle 


Table  1.  Incompressible  Parameter  Variations  to  Assess  the  Effects  on  Velocity 


Table  2.  Compressible  Parameter  Variations  to  Assess  the  Effects  on  Velocity 
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were  analyzed.  Changes  were  made  to  these  parameters  and  the  corresponding 
peak  velocities  were  plotted.  When  examining  the  effect  of  the  parameters  only 
flow  between  the  teeth  was  plotted  because  it  represents  the  most  significant  flow 
of  fluid  and  the  highest  rates.  In  the  compressible  flow  results,  the  position  at 
which  the  flow  reached  Mach  1  was  noted  with  respect  to  changes  in  the  design 
parameters.  Also,  for  compressible  flow,  examination  of  the  temperature  and 
pressure  rises  with  percentage  of  mesh  angle  are  shown  to  give  insight  to  the 
thermal  conditions  in  the  mesh  region. 

4.2.  Incompressible  Flow 

An  example  of  the  fluid  velocities  attainable  for  the  lubricant  using  incom¬ 
pressible  flow  are  shown  in  Figure  4.1.  Tie  mesh  simulation  was  done  with  an 
identical  gear  and  pinion  with  31  teeth  and  a  diametral  pitch  of  8.5  teeth  per  inch, 
a  pitch  line  speed  of  150  ft/s,  a  pressure  angle  of  22  degrees,  suid  a  backlash  of 
.002  inches.  The  three  curves  represent  flow  out  the  ends,  between  the  teeth,  and 
combined  flow.  The  highest  rate  occurs  for  the  teeth  flow,  v  =  30  corresponding 
to  a  velocity  of  4500  ft/s,  followed  by  combined,  V  =  13  corresponding  to  1950 
ft/s,  and  end  flow,  v  —  5  corresponding  to  750  ft/s.  The  peak  in  velocity 
is  reached  at  48%  of  the  total  mesh  cycle  for  all  three  flows.  The  flow  slows 
to  almost  zero  at  70%  of  the  mesh.  The  velocity  is  represented  by  equation  2.7, 
V / rw  =  1 ! A  X  dY j d9,  where  A  is  the  discharge  area  and  dV J dO,\s  the  change 
in  volume  with  respect  to  position.  From  the  equation,  it  is  possible  to  see  that 
smaller  discharge  areas  and  higher  volume  changes  give  larger  velocities.  Teeth 
flow  gives  the  highest  velocity  rates  because  it  presents  the  smallest  discharge 
area.  The  combined  area,  equation  2.9,  is  the  next  smallest  area,  then  finally  the 
end  area,  equation  2.3,  is  the  largest  area  of  the  three.  The  likely  cause  of  the 
peak  occurring  at  48%  mesh  is  the  area  of  discharge  has  reached  its  minimum 


Figure  4.1  Dimensionless  velocity  for  the  three  flow  areas 
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valae  while  the  rate  of  volume  change  continues  to  decrease.  The  flow  slows  to 
zero  nezir  70%  mesh  because  the  change  in  volume,  dV ] d6,  has  become  almost 
negligible. 

The  influence  of  varying  the  backlash  on  peak  velocity  between  the  teeth 
values  is  shown  in  Figure  4.2.  The  backlash  was  varied  from  .0005  inches  to 
.004  inches  using  the  same  gear  and  pinion  at  a  pitch  line  speed  of  150  ft/s. 
The  backlash  value  was  added  to  the  chord  length  calculated  between  the  two 
meshing  teeth.  The  velocity  increases  with  a  decreaise  in  backlash.  Similar  results 
are  obtained  for  other  gear  geometries.  By  assuming  no  expansion  of  the  gear 
teeth  due  to  local  temperature  rises,  the  flow  discharge  area  between  the  teeth 
can  be  controlled  by  varying  the  backlash.  Referring  to  the  velocity  equation, 
equation  2.7,  it  can  be  seen  that  as  the  discharge  area  decreases,  the  inverse  of 
the  area  increases  the  velocity.  If  no  backlaish  existed,  the  discharge  area  would 
become  very  small  as  the  teeth  approach  contact  giving  infinitely  large  velocity 
before  it  finally  closed  and  stopped  the  flow. 

A  linear  curve  results  for  changes  in  pitch  line  velocity.  A  curve  of  peak 
velocity  values  for  pitch  line  speeds  ranging  from  100  ft/s  to  400  ft/s  is  shown 
in  Figure  4.3.  Even  at  low  pitch  line  speeds  of  100  ft/s,  incompressible  theory 
predicts  high  velocities  for  this  geometry.  The  curve  is  linear  because  in  the 
dimensionless  velocity  equation,  equation  2.7,  the  velocity  is  only  dependant  on 
current  geometry  values  therefore  changes  in  pitch  line  speed  do  not  affect  the 
dimensionless  velocity.  Similar  results  are  predicted  for  other  geometries. 

The  diametral  pitch  is  the  number  of  teeth  per  inch  along  the  pitch  circle. 
Pitches  less  than  16  teeth/inch  are  designated  coarse  pitch  and  those  greater 
than  20  teeth/inch  are  fine  pitch.  The  effect  of  varying  diametral  pitch  from  8.5 
teeth/inch  to  24  teeth/inch  on  the  peak  velocity  of  discharged  air  is  shown  in 
Figure  4.4.  As  the  diametral  pitch  changes  from  coarse  pitch  to  fine  pitch,  the 


Dimensionless  velocity  as  a  function  of  backlash  changes 


Incompressible  velocity  vs.  pitch  line  velocity. 
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velocity  decreases.  This  is  a  result  of  the  volume  change,  dV /dd^  being  greater 
for  coarse  pitch  gears.  The  volume  change  depends  on  the  size  of  the  meshing 
tooth  entering  the  space  between  the  teeth.  The  size  of  the  teeth  depends  on 
the  inverse  of  the  diametral  pitch.  The  addendum  is  1/P  and  the  dedendum 
is  1.25/P  so  that  a  smaller  diametral  pitch  results  in  larger  teeth.  Changing 
the  diametral  pitch  and  holding  the  number  of  teeth  the  same  does  not  change 
the  total  mesh  angle  but  does  change  the  diameter  of  the  gears.  Therefore  over 
the  same  range  of  motion  a  larger  tooth  causes  greater  changes  in  the  volume 
between  the  teeth  for  each  increment  of  motion. 

By  changing  the  drive  ratio  between  the  gear  and  pinion  the  effect  that  the 
mating  geeir  size  has  on  the  fluid  velocity  can  be  seen  in  Figure  4.5.  The  driving 
gear  had  31  teeth  with  a  diametral  pitch  of  8.5  teeth/inch,  a  backlash  of  .002 
inches  at  a  pitch  line  speed  of  150  ft/s.  The  smaller  the  drive  ratio  the  higher 
the  air  velocity.  As  the  drive  ratio  increases, (  increase  the  number  of  mating 
gear  teeth),  the  total  angle  of  mesh  increases.  Since  the  size  of  the  teeth  remain 
constant,  for  a  larger  drive  ratio  the  meshing  tooth  causes  smaller  changes  in 
volume  as  it  passes  through  the  longer  cycle. 

The  peak  velocities  for  pressure  angles  20,  22,  and  25  degrees  axe  given  in 
Figure  4.6.  Little  variation  in  peak  velocity  occurs  for  the  three  values.  This 
is  primarily  due  to  the  small  change  in  tooth  size.  Changing  the  pressure  angle 
makes  the  base  circle  smaller  for  larger  pressure  angles.  The  base  circle  is  given 
as  =  ^pitch  ^  COS(f>  where  <j)  is  the  pressure  angle. 

Changing  the  percentage  of  flow  distributed  between  the  teeth  and  ends  is 
given  in  Figure  4.7.  This  change  increased  the  discharge  area  which  slowed  the 
velocity. 


Figure  4.5  The 


ity  changes  due  to  pressure  angle  variation 
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4.3.  Compressible  Flow 

Taking  into  account  the  compressibility  of  air,  the  meshing  action  of  the 
gears  compressed  and  pumped  out  the  air  trapped  between  the  meshing  teeth. 
An  example  of  the  velocities  for  the  three  flow  areas  is  given  in  Figure  4.8.  An 
identical  gear  and  pinion  with  31  teeth,  8.5  teeth/inch  a  backlash  of  .002  inches 
were  meshed  at  a  pitch  line  velocity  of  150  ft/s.  The  rate  of  the  three  flows 
gradually  increased  from  zero  velocity.  The  position  at  which  sonic  velocities 
are  reached  can  be  seen  by  the  abrupt  change  in  the  velocity  curve.  The  flow 
between  the  teeth  was  the  first  to  reach  Mach  1  at  35%  of  the  mesh  cycle.  The 
combined  flow  eventually  reached  Mach  1  at  47%  of  the  mesh  cycle  and  the 
endflow  did  not  reach  Mach  1.  The  flow  between  the  teeth  gives  the  highest 
velocity,  v  =  8.3  which  corresponds  to  1245  ft/s,  followed  by  combined  flow, 
V  =  8  which  corresponds  to  1200  ft/s,  and  end  flow,  V  —  6  which  corresponds 
to  900  ft/s.  As  the  mesh  begins,  the  pressure  change  inside  the  tooth  space 
increases.  Once  Max:h  1  is  reached,  the  flow  is  choked  and  remains  at  Mach  1. 
Any  changes  in  velocity  are  due  to  a  local  temperature  rise.  The  difference  in 
flow  velocity  between  the  three  assumed  areas  can  be  attributed  to  the  difference 
in  discharge  areas.  A  wide  set  of  gears  would  have  the  highest  velocities  because 
the  flow  tends  to  move  out  between  the  teeth  instead  of  the  ends. 

For  the  same  set  of  gears  and  a  pitch  line  speed  of  150  ft/s,  the  effect  of 
changing  the  backlash  from  .0005  inches  to  .004  inches  is  shown  in  Figure  4.9. 
Figure  4.10  gives  the  percent  of  mesh  when  the  air  velocity  reaches  sonic  levels. 
As  expected,  the  smaller  backlash  allows  the  velocity  to  reach  sonic  levels  sooner. 
Changing  the  driven  gear  to  63  teeth  and  keeping  the  pitch  line  velocity  at 
150  ft/s  results  in  subsonic  flow  throughout  the  mesh  cycle.  Figure  4.11  gives 
subsonic  results  for  changing  backlash.  The  effects  for  compressible  flow  are  not 
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Figure  4.9  Peak  velocity  for  backlash  variation. 
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Figure  4.10  Position  Mach  1  is  reached  for  backlash  changes. 
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as  dramatic  as  the  incompressible  results,  Figure  4,2,  this  is  due  to  one  part  of 
the  pressure  term,  P JV  dV being  more  dominant  in  equation  3.16. 

For  a  31  teeth  driver  and  a  63  teeth  driven  gear  with  a  diametral  pitch  of 
8.5  teeth/inch,  the  geometry  does  not  produce  sonic  velocity  levels  for  changes 
in  pitch  line  velocity.  The  effect  is  seen  in  Figure  4.12.  As  the  pitch  line  velocity 
increases,  the  resulting  velocity  and  begins  to  level  off  near  1100  ft/s  for  pitch 
line  velocities  greater  than  350  ft/s.  The  non  ’inear  response  of  the  peak  velocity 
is  due  to  the  non-linear  pressure  rise  and  the  7  power  in  the  velocity  equation. 
A  sonic  case  is  given  in  Figure  4.13,  the  corresponding  position  at  which  sonic 
was  reached  is  in  Figure  4.14.  The  pressure  rise  is  seen  in  Figure  4.15.  From  the 
pressure  curve,  it  can  be  seen  that  the  velocity  follows  the  same  form. 

The  effect  of  varying  the  drive  ratio  by  changing  the  number  of  teeth  in  the 
driven  gear  is  shown  in  Figure  4.16.  For  large  drive  ratios,  the  air  velocity  is 
slower.  The  effect  of  drive  ratio  on  the  position  [%  mesh)  when  the  air  reaches 
sonic  is  given  in  Figure  4.17.  The  effect  is  primarily  due  to  the  longer  approach 
associated  with  a  larger  drive  ratio.  The  longer  approach  reduces  the  volume 
changes  with  respect  to  position. 

Fixing  the  number  of  teeth  at  31  and  changing  the  diametral  pitch  from 
8.5  teeth/inch  to  24  teeth/inch  gave  the  velocities  shown  in  Figiire  4.18.  The 
effect  of  changing  the  diametral  pitch  on  position  sonic  is  given  in  Figure  4.19. 
The  smaller  the  diametral  pitch  the  sooner  it  reaches  sonic.  Subsonic  geometry 
results  for  changing  diametral  pitch  are  shown  in  Figure  4.20.  When  keeping  the 
number  of  teeth  constant,  the  diametral  pitch  does  not  change  the  total  angle 
of  approach,  but  does  change  the  diameter  of  the  gear.  This  means  that  for  the 
same  angle  of  approach  the  difference  in  velocities  is  primarily  due  to  the  change 
in  the  volume  with  respect  to  position,  dV J dQ .  The  major  factor  affecting  the 
change  in  volume  is  the  size  of  the  meshing  tooth.  The  tooth  size  varies  inversely 
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Figure  4.14  Mesh  position  at  which  Mach  1  is  reached  for  pitch  line  changes. 
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Figure  4.16  Peak  velocity  due  to  drive  ratio  changes  (changes  in  the 
number  of  driven  pinion  teeth  )• 
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with  the  diametral  pitch.  Smaller  diametral  pitches  result  in  larger  teeth  which 
cause  greater  changes  in  volume  for  each  incremental  step. 

The  effect  of  changing  the  pressure  angle  from  20  to  22  and  then  25  degrees 
gives  the  result  in  Figure  4.21.  A  slight  variation  occurs  due  to  the  small  change 
in  tooth  size  for  changing  pressure  angle. 

A  typical  pressure  variation  for  the  three  flows,  end,  teeth  and  combined, 
is  given  in  Figure  4.22.  The  pressure  peaks  correspond  to  the  peak  velocity 
position. 

Changing  the  flow  distribution  by  making  the  endflow  area  account  for  a 
larger  percentage  of  the  total  disch::rge  area  is  given  in  Figure  4.23.  As  expected, 
as  greater  percentages  of  the  end  area  axe  added  the  discharge  area  increases  and 
results  in  a  decrease  in  peak  velocity. 

The  temperature  changes  inside  the  meshing  region  of  the  teeth  for  the  three 
flow  areas  is  seen  in  Figure  4.24.  Starting  at  the  ambient  temperature,  80  F, 
the  end  flow  temperature  problem  rises  to  150  F,  the  teeth  flow  to  350  F,  and 
the  combined  flow  problem  to  220  F.  The  energy  generated  must  be  dissipated. 
A  portion  of  the  energy  is  dissipated  in  the  air-oil  mix  and  the  other  portion  is 
absorbed  by  the  gear  teeth.  If  too  much  energy  is  absorbed  by  the  teeth  thermal 
distortion  results. 

4.4.  Comparison  of  Incompressible  and  Compressible  Flow 

Using  identical  gear  and  pinion  of  31  teeth  and  8.5  teeth/inch  diametral  pitch, 
the  velocities  for  incompressible  and  compressible  flow  were  compared  for  each 
of  the  assumed  discharge  flow  areas  at  a  pitch  line  speed  of  150  ft/s. 

Comparing  the  end  flow  in  Figure  4.25,  the  incompressible  reaches  its  peak 
at  48%  of  the  mesh  cycle.  The  compressible  flow  case  reaches  a  peak  at 
approximately  60  %  of  the  mesh  cycle  while  both  have  approximately  the  same 


Dimensionless  pressure  in  the  mesh  region  for  the  three  flow  areas, 


4,24  Temperature  in  Fahre^.r:  (  ir  the  mesh  region  for  the  three  flow  areas 
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peak  value.  The  delay  in  the  compressible  peak  occurs  because  of  the  compressing 
of  the  air  and  the  gradual  rise  in  pressure.  Comparing  the  flow  between  the  teeth, 
Figure  4.26,  the  incompressible  peaks  at  V  =  30  and  the  compressible  reaches 
Mach  1  at  35  %  of  the  mesh  cycle.  A  direct  comparison  of  peak  velocities 
for  incompressible  and  compressible  flow  is  not  meaningful  because  once  the 
compressible  velocity  reaches  Mach  1  it  is  choked  and  remains  at  Mach  1. 

Comparing  the  combined  flow  in  Figure  4.27  the  incompressible  is  again 
higher  and  the  compressible  reaches  sonic  levels  at  43  degrees. 

4.5.  Conclusions 

A  one-dimensional  approximation  of  the  air  oil  flow  resulting  from  the  mesh¬ 
ing  of  spur  gear  teeth  indicates  the  velocity  reaches  high  rates.  Incompressible 
theory  predicts  oil  velocities  as  high  as  4500  ft/s  depending  on  the  pitch  line 
speed.  Using  compressible  flow  theory,  some  geometries  produce  air  velocities 
high  enough  to  approach  and  reach  the  speed  of  sound,  Mach  1,  at  the  exit 
plane.  If  the  velocity  reaches  Mach  1,  the  conditions  may  lead  to  the  formation 
of  a  shock  wave  downstream  of  the  exit  plane.  Should  a  shock  wave  form,  noise 
is  emanated  which  contributes  to  the  overall  noise  of  the  gear  mesh. 

Since  the  velocities  were  shown  to  depend  on  a  number  of  different  design 
parameters — diametral  pitch,  backlzish,  pitch  line  velocity,  pressure  angle,  drive 
ratio,  and  face  width — further  analysis  should  be  conducted  to  quantify  their 
exact  effect.  This  thesis  has  shown  the  parameters  most  significantly  affecting 
the  air  velocity  appear  to  be  the  drive  ratio  and  the  face  width. 

Further  considerations  of  this  problem  should  include  an  experimental 
analysis  of  the  fluid  flow.  Also,  an  expansion  of  the  analytical  portion  of  the 
problem  should  include  two  and  three  dimensions  to  determine  the  flow  pattern 
for  different  width  gears  The  inclusion  of  the  effects  two  rotating  teeth  have  on 
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the  velocity,  the  effect  of  the  moving  exit  plane,  and  a  mixture  of  air-oil  mist 
velocity  would  make  the  problem  closer  to  the  actual  system. 


Figure  4.26  Comparison  of  incompressible  and  compressible  teeth  flow 
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Figure  4.27  Comparison  of  incompressible  and  compressible  combined  flow. 
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APPENDIX 

COMPUTER  ALGORITHM  GUIDE  AND  A  FLOW  CHART 

1.  PROGRAM  EXPLANATION 

The  computer  algorithm  approximates  the  velocity  for  flow  between  meshing 
gear  teeth  by  assuming  incompressible  and  compressible  flow  theories.  For  each 
case,  three  flow  areas  were  used:  end  flow,  teeth  flow,  and  combined  flow.  The 
program  was  written  in  Fortran  '/  7  on  the  Data  General  Eclipse  system.  It  is  run 
interactively  by  specifying  the  data  file  and  the  output  listflle. 

2.  MAIN  PROGRAM 

The  main  program  reads  and  echoes  the  input  information.  The  driven 
teeth  were  considered  the  pinion  and  the  gear  was  the  driver.  The  first 
subroutine  called  is  SUBROUTINE  GEOCAL,  which  computes  the  gear  and 
pinion  geometry.  For  the  chosen  set  of  gear  geometries,  the  intial  contact 
(starting  position  gear  tooth  centerline)  is  computed.  The  centerline  position 
is  incremented  in  a  DO  LOOP  by  the  chosen  step  size.  The  program  stops  when 
full  mesh  position  is  reached. 

The  subroutine  to  compute  the  incompressible  velocity  is  SUBROUTINE 
INCOMPRESSIBLE.  The  value  for  the  gear  centerline  position,  its  absolute 
value,  and  the  step  size  is  transferred  to  the  subroutine.  Values  for  the  three 
flow  area  velocities  are  returned. 

For  compressible  flow  computations,  SUBROUTINE  COMPRESSIBLE  is 
called  to  compute  the  pressure,  velocity,  temperature,  and  Mach  #  for  each  of 
the  three  flow  areas.  The  position  of  the  gear  centerline  and  the  size  of  the  step 
are  sent  into  the  subroutine. 
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The  computed  values  are  written  to  the  listfile  (default  unit=12)  at  the 
desired  frequency  controlled  through  the  input  file.  The  position  at  which  Mach  1 
is  reached  is  also  written  onto  the  listfile. 

3.  SUBROUTINE  AREA 

This  subroutine  is  called  by  SUBROUTINES  INCOMPRESSIBLE  and 
COMPRESSIBLE,  it  computes  the  area  between  the  meshing  teeth,  AIR2, 
by  first  determining  the  chord  length  between  the  teeth,  CHORD2,  and  then 
multiplying  by  the  constant  face  width.  The  chord  is  computed  by  locating  the 
trailing  tip  of  the  driving  tooth  and  computing  the  arc  to  the  driven  tooth.  The 
backlash  amount  is  added  in  the  incompressible  and  compressible  subroutines. 

4.  SUBROUTINE  GEOCAL 

This  subroutine  uses  the  standard  AGMA  and  USASI  tooth  systems  for  spur 
gears.  The  computed  radii  and  tooth  widths  are  written  to  the  output  listfile. 

5.  SUBROUTINE  DERIVATIVE 

This  subroutine  computes  the  change  in  volume  with  respect  to  position, 
dV  j d9.  It  uses  the  forward  difference  method  for  the  derivative  at  the  starting 
position  and  the  remaining  derivatives  are  computed  by  the  central  difference 
method.  It  calls  the  SUBROUTINE  VOLUME. 

6.  SUBROUTINE  COMPRESSIBLE 

This  subroutine  computes  the  pressure,  temperature,  velocity,  and  mach 
number  for  the  three  flow  areas  by  assuming  compressible  flow.  It  is  appropriate 
for  lightly  lubricated  gears  where  the  primary  discharged  fluid  is  air.  It  uses  the 
Runge-Kutta  fourth  order  formula  to  solve  the  change  in  pressure  with  respect 
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to  position,  dP  j dQ.  The  subroutine  calls  SUBROUTINE  RHS  to  compute  the 
intermediate  values  of  the  function  for  the  Runge-Kutta  method.  The  formula  is 

Vj+v  =  Vi  + 

where 

Ac 

=  Vj  + 

yj+^  ~  yj 

2/y+i  =  yy  +  ^^/(yy+^»*y+^) 

The  intermediate  values  2/v  ,  i,  2/vT.i.  j  o,Tid  J/v  .  ,  must  be  computed  in  the  order 
given  since  they  are  interdependent  [15], 

7.  SUBROUTINE  RHS 

This  subroutine  computes  the  right  hand  side  of  the  pressure  equation 
for  the  three  flow  areas.  It  calls  SUBROUTINES  AREA,  VOLUME,  AND 
DERIVATIVE. 

8.  SUBROUTINE  INCOMPRESSIBLE 

This  subroutine  computes  the  velocity  for  the  three  flow  areas  by  assuming 
an  incompressible  fluid.  This  case  is  appropriate  for  heavily  lubricated  gears. 
SUBROUTINES  AREA  and  VOLUME  are  called. 

9.  SUBROUTINE  BISEC 

This  subroutine  determines  the  root  of  an  equation  that  falls  within  the  given 
limits.  It  is  called  by  the  SUBROUTINE  VOLUME  to  determine  the  point  where 
the  involute  portion  of  the  driving  gear  intersects  the  addendum  circle  of  the 


pimon. 
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10.  SUBROUTINE  VOLUME 

This  subroutine  computes  the  volume  of  space  between  the  driver  and  the 
driven  teeth.  The  subroutine  is  broken  up  into  various  integral  sets  to  compute 
the  area  of  the  driving  geaj  that  overlaps  into  the  meshing  region.  This  area 
is  used  to  determine  the  two-dimensional  plane  area  which  is  multiplied  by  the 
constant  face  width  to  give  a  volume.  SUBROUTINES  BISEC  and  SUBTRACT 
are  called. 

11.  SUBROUTINE  APPROX 

This  subroutine  determines  the  end  and  beginning  points  for  the  cubic 
fit  to  the  approximation  for  the  area  to  subtract  when  the  SUBROUTINE 
SUBTRACT  method  is  no  longer  valid.  The  cubic  fit  used  was  f{x)  = 
Co  +  aix  -f  ajX'  +  a^x^. 

12.  SUBROUTINE  SUBTRACT 

This  subroutine  computes  the  amount  of  area  to  be  excluded  from  the  2-D 
plane  area  when  determining  the  volume. 

13.  LISTFILE 

The  input  data  is  echoed  first  followed  by  the  computed  gear  and  pinion 
geometry.  The  remainder  of  the  output  consists  of  the  computed  velocities. 
The  first  column  is  the  percent  mesh  position.  The  incompressible  flow  columns 
are  designated  INC  E,  INC  T,  INC  C  and  the  compressible  flow  columns  are 
COMP  E,  COMP  T,  COMP  C  and  MACH  E,  MACH  T,  MACH  C  where  E 
corresponds  to  end  flow,  T  for  teeth  flow,  and  C  for  combined  flow. 
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14.  INPUT  FILE 

The  following  variables  are  read  free  format  one  per  line  into  the  program: 
number  of  gear  teeth 
pitch  of  the  gear 

pressure  angle  of  the  gear  in  degrees 

pressure  angle  of  the  pinion  in  degrees 

pitch  of  the  pinion 

number  of  pinion  teeth 

coarse  or  fine  pitch  gear  (coarse=l,fine=2) 

coarse  or  fine  pitch  pinion  (coarse=l,fine=2) 

the  maximum  iterations  in  the  bisec  routine 

the  convergence  criterion  for  bisection 

face  width  of  gear  in  inches 

face  width  of  pinion  in  inches 

convergence  criterion  for  derivative  calculations 

density  of  air  in  /fern/ 

pressure  of  the  ambient  Ibf  /iri^ 

dimensional  gravitational  constant  in  Ib^injlb 

pitch  line  velocity  ftj S 

ratio  of  specific  heats 

increment  of  change  A0  in  degrees 

total  number  of  loops  executed 

counter  for  printing  results 

the  number  ot  iterations  between  write  statements 


backlash  in  inches 
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